Perturbative Renormalization Factors of Baryon Number 
Violating Operators for Improved Quark and Gauge Actions in 

Lattice QCD 



Sinya Aoki^, Yoshinobu Kuramashi*^ *, Tetsuya Onogi'^ and Naoto Tsutsui'^ 
""Institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305-8571, Japan 
''Department of Physics, Washington University, St. Louis, Missouri 63130, USA 
"^Department of Physics, Hiroshima University, Higashi- Hiroshima, Hiroshima 739-8526, Japan 

(February 1, 2008) 

Abstract 

We calculate one-loop renormalization factors of three-quark operators, 
which appear in the low energy effective Lagrangian of the nucleon decay, 
for 0(a)-improved quark action and gauge action including six-link loops. 
This calculation is required to predict the hadronic nucleon decay matrix ele- 
ments in the continuum regularization scheme using lattice QCD. We present 
detailed numerical results of the one-loop coefficients for general values of 
the clover coefficients employing the several improved gauge actions in the 
Symanzik approach and in the Wilson's renormalization group approach. The 
magnitudes of the one-loop coefficients for the improved gauge actions show 

sizable reduction compared to those for the plaquette action. 
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I. INTRODUCTION 



While nucleoli decay is one of the most exciting prediction from grand unified theo- 
ries (GUTs) with and without supersymmetry, none of the decay modes have been exper- 
imentally detected up to now. Furthermore, the ongoing Super-Kamiokande experiment is 
now pushing the lower limit on the partial lifetimes of the nucleon by an order of magni- 
tude from the previous measurements. In principle this would give a strong constraint on 
(SUSY-) GUTs, however, the uncertainties in the theoretical prediction of lifetimes due to 
poor knowledge of quantum effects at low energy such as hadron or SUSY scales obscure 
the direct impact of the experimental lifetime bound on the physics at the GUT scale. In 
particular, one of the main uncertainties has been found in the evaluation of the hadron 
matrix elements for the nucleon decays, for which various QCD models have given estimate 
differing by a factor of ten. Therefore precise determination of the nucleon decay matrix 
elements from first principles is required, for which Lattice QCD can play a crucial role. 

Recently we carried out a model-independent calculation of the nucleon decay matrix 
elements employing the Wilson quark action and the plaquette gauge action in the quenched 
approximation [jl| . Although one naively expects the error in the discretization and quench- 
ing approximation which are the two main systematic errors, it would be desirable to reduce 
these unknown systematic errors for high precision calculation. As a step toward this goal 
we are required to reduce the scaling violation effects by improving quark and gauge actions. 
In this article we present perturbative results for renormalization factors of baryon number 
violating operators for improved quark and gauge actions: the 0(a)-improved "clover" ac- 
tion originally proposed by Sheikholeslami and Wohlert and the gauge action improved 
by addition of six link loops to the plaquette term in the Symanzik approach and in 
the Wilson's renormalization group approach 0,0. Values of the one- loop coefficients of 
the renormalization factors are numerically evaluated for combinations of general values of 
the clover coefficients in the quark action and some specific values of the coefficients of the 
six-link loop terms in the gauge action. 
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This paper is organized as follows. In Sec.|| we give the improved quark and gauge 
actions on the lattice and their Feynman rules relevant for their calculation. Our calcula- 
tional procedure of the renormalization factors for the baryon number violating operators 
is described in SecPTT], where we present expressions and numerical values for the one- loop 
coefficients of the renormalization factors. Our conclusions are summarized in Sec.|V|. 

The physical quantities are expressed in lattice units, and the lattice spacing a is sup- 
pressed unless necessary. Throughout this paper we use the same notation for quantities 
defined on the lattice and their counterparts in the continuum. In case of any possibility of 
confusion, however, we shall make a clear distinction between them. 



II. ACTIONS AND FEYNMAN RULES 



For the gauge action we consider the following general form including the standard 
plaquette term and six-link loop terms: 
1 



•S'gauge = "( Co '^^^P' + '^1 Tr?7rtg + C2 ^ TrUchr + C3 

3 I nlannftttp rftrtariErlfi chair 



pig 



plaquette rectangle 

with the normalization condition 



parallelogram 



Co + 8Ci + 16C2 + 8C3 = 1, 



1) 



(2) 



where six-link loops are 1x2 rectangle, a bent 1x2 rectangle (chair) and a three-dimensional 
parallelogram obtained by multiplying the link variables 



Un,, = exp (^tg Y T^A^in + f,/2)j . 



The free gluon propagator is obtained in Ref. [^: 
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2\2 



with 



k^ = 2sin f Y 



(3) 

(4) 

(5) 
(6) 



The matrix Afj_„ satisfies 



(i) A^^ = for all (7) 

(ii) A^^ = A^^, (8) 

(iii) VW= (9) 

(iv) V(0) = 1 for /X ^ z/, (10) 



and its expression is given by 



ZA4 



-\-(k (?l/p?;/T^;/ + QvpQpTkp + Q'i/-7-?pT^r) 

+%p(lur{kl + kl){kl + kl) + g^.g.p(A;2 + + 
~Q.tiv(lpT{k^p + ^r)^ ~ {QupQvp + (liJ.TQ.vT)kpk^ 

~Qtj.v{QtJ.pk^k'^ + Qprk^k'^ + Qupklk'^ + Qurklk'^) , (11) 



with II ^ V ^ p ^ T the Lorentz indices. 5^1^ and A4 are written as 



= (1 - 5^,) [1 - (Ci - C2 - C^){kl + ^^) - (C2 + c^)P\ , (12) 



^4 = 13 ]J + X] klklqi,^{q^pq^r + gMr^i/p)- (13) 

In the case of the standard plaquette action, the matrix ^4^^ is simplified as 

^plaquette ^ ;l - 6^,. (14) 

For the quark action we consider the 0(a)-improved quark action: 

-^quark = J2 nl2 {^n{-r + -fi,)Un,p,llJn+fl + i^n{-r - Tm)^^!-^,//'^^-^} + ("^0 + 4r) ^^ni^n 
n fi n 

-CswY.Y.^dji'nCTnuPp.Mlpn, (15) 
n ix,u ^ 

where we define the Euclidean gamma matrices in terms of the Minkowski matrices in the 
Bjorken-Drell convention: 7^- = -i'j^BD U = 1; 2, 3), 74 = 7^^, 75 = Ybd and a^,, = ^[7^, 'f^]. 
The field strength P^^^ in the "clover" term is given by 



Ui{n 
U2{n 

UJn 



TT TT - TT^^ TT'^ 

^n,fj.'-^ n+fi,u'-^ n+u,ti^ n.u^ 
U U'^ . . . TT - 

'-^n-f^,^l'-^n-^J,-l/,I/'-^rl-fl-u,fl'-^n-l/,l/^ 
TT'^ TT - TT ~ ~ TT^ 



From the quark action (p!5D we obtain the free quark propagator 
Sq^{p) = iX!7A.sin(A;^) +mo + r^(l -cos(p^)). 



(16) 
(17) 
(18) 
(19) 
(20) 

(21) 



In order to calculate renormalization factors of the baryon number violating operators up 
to one-loop level, we need the following vertices, 



(l) = -gT^Cs^^- ( ^/«^sin(p^ - q^) J cos ^ ^ 



— — zrsm — 

2 J V 2 



(22) 
(23) 



with incoming quark momentum and outgoing quark momentum. The first vertex 
originates from the Wilson quark action and the second one is the interaction due to the 
clover term. In the present calculation of the vertex corrections for the baryon number 
violating operators, the two-gluon vertices with quarks give no contribution. 

We should note that the baryon number violating operators contain a charge conjugated 
field, whose action is obtained from eq.(|15|) with the replacement of 



-{T 



a\T 



(24) 



where the superscript T means the transposed matrix. According to this change, the Feyn- 
man rule of the quark-gluon vertices in eqs.(p2D and (^) should be modified for the charge 
conjugated field. 



III. RENORMALIZATION FACTORS FOR BARYON NUMBER VIOLATING 

OPERATORS 
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A. Calculational procedure 



We consider the following baryon number violating operators in the continuum and on 
the lattice: 

{OZy), = ^""^ [{riTTx{iJ2?] [ry(^3)l5 , (25) 
{O'k'k), = e""'^ [{1 + rm,{l - z)} {r,TVx{i^2f 



-z'-{riprvx{m' 



X 



1 + ^mo(l - ^)| Vyi^^y - 4ry(f03) 



(26) 



with 



<J = ^ E - <-a^:-A,m) 7/., (28) 



where il^'^ = C with C = 7472 is a charge conjugated field of ip. Dirac structures are 
represented by Tx <8) Fy = Pr Pr, Pr (g) P^,, Pl <S> Pr, Pl ® Pl with the right- and left- 
handed projection operators Pr^l = (1 =t 75)/2. The summation over repeated color indices 
a, b, c is assumed. 

Ultraviolet divergences of composite operators are regularized by the cutoff in the 
lattice regularization scheme, while this is achieved by a reduction of the space-time di- 
mension from four in some continuum regularization schemes, where we consider the naive 
dimensional regularization (NDR) scheme and the dimensional reduction (DRED) scheme. 
Operators defined in different regularization schemes can be related by renormalization fac- 
tors: 

O'^^fil,) = Zdiag(/^a)O^^U«) + Z^^&^^ia) (29) 
with jj, the continuum renormalization scale. The explicit chiral symmetry breaking due to 



the Wilson term in the quark action (|T5]) causes the mixing between operators with different 
chiral structures, which is denoted by &^y. 

Since QCD is a asymptotically free theory, Zdiag.mix are expected to be perturbatively 
calculable in terms of the coupling constant j (IGvr^) at high energy scales, which gives the 
following expressions, 

^3 / . „ 8 



Zdmg^l^o) = 1 + 2 \ ~ ~ 3 + ^i^J + 8 ln(/ia) + Ay^diag 



(30) 



where Sq denotes the additive mass renormalization on the lattice, is a contribution from 
the wavefunction and Ay is from the vertex function. Eq and are obtained by calculating 
the continuum and lattice quark self-energies and Ay is from the continuum and 

lattice vertex functions A™°*'^'^". 

The quark self-energies in the continuum and on the lattice are defined through the 
inverse full quark propagators: 

{Sin"\p) = + m - Yf^S™'^*(p), (32) 
{S':'r\p) = ^E7.sin(p,) + mo + r $:(! - cos(p,)) - rf^S'^^^p), (33) 

where we consider massless quark. The difference between and S^****, which originates 
from the regularization scheme dependence of the self-energy, gives the quark wavefunction 
renormalization factor, 

9Si^"(p) 



8, , , , dE^^'^Hp) 
- - ln(/ia) + A^ ^"^^ 



3 idp 
The additive quark mass renormalization is expressed as 



(34) 



p=0 



"^0 - 77-^^0 (35) 

with 

Eo = S^^"(p = 0). (36) 



Notice that the renormahzation factor (^) is given for the case of mo = (7^/(167r^)So, where 
we consider the renormahzation for massless quark. 

The vertex functions up to one-loop level in the continuum and on the lattice are ex- 
pressed in the following way, 

^cont ^ ^abc^^ ^ ^ tI^C^'Tx ® Tyl/^^, (37) 



167r 

2 

2 T A"latt 
XY 



Ibvr^ L \ ' 

+fx ® fy {y^£^ + zV^f - ^'V^^lf )] , (38) 



where the number of prime in the superscript of the lattice vertex corrections denotes the 
number of covariant derivative applied to the quark fields at the vertex. Tx ® Ty term 
represents the mixing contribution. The difference between A™y and A^*y leads to 



8 ln(/ia) + Ay,diag = V^y - (^^y + zV^f - z Vj^|f ) , (39) 
Ay,^i. = - [V'^^ + zV^f - z^Vi^'p') . (40) 

We note that the lattice quark-self energy and the lattice vertex corrections are general 
function of the clover coefficient csw in the quark action and the six-link loop parameters 
Ci 2,3 in the gauge action. Calculation of was already carried out in Ref. employing the 
general values for csw- For ci,2,3 they choose some specific values: ci = — 1/12,C2 = C3 = 
in the tree-level Symanzik improvement, ci = —0.252, C2 = 0,C3 = —0.17 suggested by 
Wilson based on renormahzation group improvement and Ci = — 0.331, C2 = C3 = and 
ci = —0.27, C2 + C3 = —0.04 by Iwasaki. According to Ref. p we evaluate Ay for general 
values of csw and for the specific values of Ci,2,3 that they employed. 

B. Vertex corrections 

We calculate the vertex corrections of the operators in eqs.(PBD and in the Feynman 
gauge employing the massless quarks and and the massless charge conjugated quark with 
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momenta pi = P2 = = as external states. The infrared singularities are regularized by 
a fictitious gluon mass A introduced in the gluon propagator. 

One-loop vertex corrections on the lattice are illustrated in Fig.^ We find that the 
lattice vertex corrections is a second polynomial function of the clover coefficients csw- The 
relevant diagrams for Vxy Vxy ci-re Figs.|I|(a) — (i), the sum of which gives 



Tfi/L ® TyVr/l,Y + ^R/L ® TyVr/l,Y 



i=0.1,2 





1 


(C^61n 


A2a2 



diag 



+ ( Tl/R ® Ty ± - ^ 7^75 (g) rY7^ I J2 Csw 



(41) 



M / j=0,l,2 

with = (A^+ 1)/ (2A^) in the SU(A^) group. The explicit forms of v^^^^ and f^j^ ^^re given 
by 



JO) 

diag 



,(1) 
diag 



,(2) 
diag 



,(0) 



,(1) 



,(2) 



r ^ {48(A3 + r'Air + 84 + 16/,} - Bin' - 



6 



(k 



2\2 



'Cjf 6 In 



TT 



- d^k 

-TT 7r2 

- d^k 

- TT TT 

- d^k 

^2 

- d^k 

-TT TT 

-TT TT 



-^{l6r^A?(A3-4A;',o)} 



(42) 
(43) 

(44) 

(45) 

(46) 

(47) 



where 



Go = 

4 = A^^i - 4A^ + (I6A3 - 4sin2(A;^))A^,o, 
h = -Al, + AI + 4(-A3 + sin2(A;^)) A^'^o, 



(48) 
(49) 
(50) 
(51) 
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^1,1 = El^A'- + A^,)sm\k^)sm\K), 

V 



(52) 
(53) 
(54) 

(55) 



We do not take the sum over the index for Aj'q- In a similar way we obtain the 

expressions of Vxy and Vxy from Figs.|I](a) — (i): 



where 



R/L 



diag 



1=0,1,2 



+ I Tl/R ® Ty ± - ^ 7^75 ® ry7^ ) ^ Csw 



i=0,l,2 



„/(0) 
^diag 



r d^k 1 r 9 . 9.1 



{l2r2Ai(A3-4A^,o)} 



FqGq 



7T^ 



diag 



V. 



V 



/(2) 
diag 

'(0) 



+ 



IT Tt 



\- {-8(r^A3 - r^A?)/4 



- d^fc r 1 



7^2 
2" 



{3r^/a} 



^ {-Sr^A.As/.} 



FqGo 



/(I) 



,'(2) 



^2 

-TT TT 



{8r^Ai(A3-4A^,o)} 



^{-8(.^A3-.^AD/4 + ^{2r^/4 



(56) 



(57) 



(58) 
(59) 



with Cf = {N^ - 1)/{2N) in the SU(A^) group. The sum of Figs.|(a)-(k), which 
the tadpole diagrams at the vertex, yields Vx y^ and V^''y*, 



(60) 
(61) 
(62) 
include 
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R/L 



Y 



^diag 



1=0.1,2 



^L/R ® Ty ± - ^ 7^75 (g) Ty^f^ J2 ^SW 



(63) 



i=0,l,2 



where 



„/'(o) 

^diag 
//(I) 



V 



V 



//(2) 
diag 

m 



+ 



1 



FqGq 



,"(1) 



,"(2) 



^ "^'^ ^ It {8r^A,A34} + -i^ {-2r^A,4} 
{4r^A?A34} 



(64) 
(65) 
(66) 



FqGq 



(67) 
(68) 
(69) 



(i) 

We present numerical values of f diag mix Table 



,g,mix 



V. 



'(i) 

diag, mix 



in Table |I| and WdtLmix in 



Table |T|, which are evaluated with r = 1 using the Monte Carlo integration routine BASES 
for specific values of Ci and C2 + C3. The numerical accuracy is better than 0.01%. 
Numerical values for fdiig,mix5 "^cUag.mix ^-^d fdiag,mix can be also obtained by using the results 
for vertex corrections of bilinear operators in Ref. 0, in which the numerical values are 
evaluated in a different way. This is used as a check of our calculation. We note that a 
special case of csw = and Ci = C2 + C3 = represents combination of the Wilson quark 
action and the plaquette gauge action, for which perturbative renormalization factors for 
the baryon number violating operators has been already calculated H,|l|]. Comparison of our 
results with theirs gives us another check of our calculation. Numerical values in Tables |, H 



and |ITl| show that the one-loop coefficients in the vertex corrections diminishes by 10 — 20% 
for the tree-level Symanzik action compared to those for the plaquette action. Further 
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reduction of the magnitude is observed for the renormahzation group improved actions. 
These features are also found in the case of bihnear operators 0. 

In the continuum, the vertex correction at one-loop level is expressed as 

,2 



Vr/lx = 6 In 
The finite constant fdiag is given by 



A2 



+ ^^diag- (70) 



<^ = I (71) 



3 

= 4 (72) 



for the NDR and DRED schemes with MS subtraction. 

The vertex corrections on the lattice give the expression of O^xy 6q.(|2 

1 



/s^latt ( \ abc 



{riTTL,R{^2?\ [ry (V^3)1 ± T E '^"^ [(^i')"7^75(V^2)'J [TyiM^Y] ■ (73) 



Comparing the results for the vertex corrections in the continuum and on the lattice, we 
obtain the vertex correction components in the renormahzation factors of eqs.(|SUD and (Plj), 

. NDR.DRED _ NDR.DRED _ „ , 'W _ z^./'^^ (7A\ 

^y.diag — ^diag Csw l^diag + ^^diag ^^diagjj \'^) 

j=0,l,2 

^VM^ = - E csw(^n^S. + ^41-^^^:22), (75) 

i=0,l,2 

where Ay mix is independent of the renormahzation scheme in the continuum. To obtain 
the diagonal part of the renormahzation factor in eg. (|30|) , we also need the wavefunction 
component A^. This quantity is already evaluated in Ref. |^ employing the NDR scheme 
in the continuum, where Cpz^ in their notation corresponds to our A^°^. We note that 
A^^^° can be obtained from Aj^°^ by 



^DRED ^ ^NDR _ 4 ^^g^ 



For the mixing part of the renormahzation factor in eg . (^I]) , A^ has no contribution. 
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IV. CONCLUSIONS 



In this paper we have calculated the one-loop contributions for the renormalization fac- 
tors of the three-quark operators employing the improved quark and gauge actions. Detailed 
numerical values of the one-loop coefficients are presented for general values of the clover 
coefficients and for some specific values of Ci, C2 and C3 in the gauge action improved by 
the Symanzik and renormalization group approaches. We find that the magnitude of the 
one-loop coefficients are considerably reduced for the improved gauge actions compared to 
those for the plaquette action, which is a desirable feature in the practical implementation 
of lattice QCD. 
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FIGURES 
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TABLES 

TABLE I. Finite parts u^iag and Vmix for vertex correction of the baryon number violating 



operators. Coefficients of the term csw" {n = 


0,1,2) are 


given in the column marked 


as (n). 


gauf 


;e action 




^diag 












Co + C"i 


fO) 


fl) 


(2) 


fo) 


fl) 


(2) 








11.084 


3.3282 


1.1498 


-3.2137 


3.3282 


-1.1386 


-1/12 





10.743 


2.9496 


0.95371 


-2.7085 


2.9496 


-1.0372 


-0.331 





10.052 


2.3000 


0.65335 


-1.9381 


2.3000 


-0.84886 


-0.27 


-0.04 


10.140 


2.3654 


0.67777 


-2.0107 


2.3654 


-0.86909 


-0.252 


-0.17 


10.020 


2.2432 


0.61679 


-1.8783 


2.2432 


-0.83258 
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TABLE II. Finite parts v^j^^g and v[^^^ for vertex correction of the baryon number violating 
operators. Coefficients of the term csw" ('^ = 0, 1, 2) are given in the column marked as (n). 



gauge action v'^^^^ v'^^ 



Ci Co + C?, 


(0) 


(1) 


(2) 


(0) 


(1) 


(2) 








-25.821 


13.753 


-1.1792 


-6.3894 


4.5991 


-1.1792 


-1/12 





-21.129 


11.779 


-1.0080 


-5.0767 


3.8148 


-1.0080 


-0.331 





-14.319 


8.6222 


-0.72746 


-3.2392 


2.6134 


-0.72746 


-0.27 


-0.04 


-14.873 


8.8984 


-0.75215 


-3.3803 


2.7111 


-0.75215 


-0.252 


-0.17 


-13.583 


8.2648 


-0.69439 


-3.0321 


2.4672 


-0.69439 
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TABLE III. Finite parts v^'j^^g and v^j^ for vertex correction of the baryon number violating 
operators. Coefficients of the term csw" {n = 0,1,2) are given in the column marked as (n). 



gauge action v'^^^^ 



Cl 


('2 + C.i 




(1) 


(2,) 




(1) 


(2) 








-1.3783 


1.1792 


-0.35763 


4.5459 


-2.6297 


0.54935 


-1/12 





-1.1210 


1.0080 


-0.31313 


3.5094 


-2.0995 


0.44985 


-0.331 





-0.74755 


0.72746 


-0.23554 


2.1314 


-1.3427 


0.30008 


-0.27 


-0.04 


-0.77933 


0.75215 


-0.24280 


2.2182 


-1.3932 


0.31101 


-0.252 


-0.17 


-0.71136 


0.69439 


-0.22656 


1.9313 


-1.2279 


0.27829 
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